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I evaluate the current results of spin-dependent weakly interacting massive particle (WIMP) 
searches within a model-independent framework, showing the most restrictive limits to date derive 
from the combination of xenon and sodium iodide experiments. The extension of this analysis to 
the case of positive signal experiments is elaborated. 
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It is well known that the spin-dependent WIMP- 
nucleus interaction is quenched by spin pairing effects 
in the nucleus. For this reason, ideal nuclei for a spin- 
dependent WIMP search would be both odd Z and 
odd neutron number (N), hence even mass number A. 
Unfortunately, actual detectors use only odd A nuclei, 
so that their response is dominated by the WIMP- 
proton (WIMP-neutron) interaction only. Interpretation 
of the results has customarily proceeded by assuming 
the WIMP interacts with only the dominant group of 
nucleons. In this "odd group" approximation, an ex- 
periment using only odd Z isotopes cannot constrain a 
theory predicting a WIMP which couples with neutrons 
only, i.e., the exclusion results quoted by such experi- 
ment are WIMP model dependent. This is particularly 
true with the current WIMP candidate: the neutralino is 
a superposition of higgsinos and gauginos, whose weight- 
ings determine the coupling strengths. 

On the other hand, some odd Z nuclei, like F and I, 
exhibit a non-negligible neutron group spin, making it 
reasonable to go beyond the "odd group" approxima- 
tion. Recently, some experiments began to adopt this 
approach 0, H 0, E3 , usm S (except for DAMA) a full 
spin-dependent treatment proposed by Tovey et. al.|j|. 

Here, the impact of this framework on the interpre- 
tation of current experiments is considered, first in the 
case of null results and then for the case of a positive 
WIMP signal detection such as DAMA/Nal. This letter 
focuses on experiments that have already published spin- 
dependent limits; others will be considered in a forth- 
coming paper. The combination of xenon- and sodium 
iodide-based experiments is seen to provide the most re- 
strictive current limits on the existence of spin-dependent 
WIMP dark matter, with the fluorine-based (F-based) 
experiments offering a good possibility for near-term im- 
provements 0. 

At tree level, the general (zero momentum transfer) 
WIMP-nuclide spin-dependent elastic scattering cross 
section a a for a nucleus of mass number A is El : 



where < S PjTl > are the expectation values of the proton 
(neutron) groups spin, Gf is the Fermi coupling con- 
stant, a p<n are the effective proton (neutron) coupling 
strengths which appear in the effective lagrangian, fiA is 
the WIMP-nuclide reduced mass, and J is the total nu- 
clear spin. A discussion on the generality of Eqn. can 
be found in Ref. 0- When limits on a a are measured, 
experiments usually quote limits only for one nucleon. 
This is done by calculating, for each sensitive nuclide, 
one of the following cross sections: 
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where tr^ m is the upper limit on a a obtained from experi- 
mental data, n Pt n is the WIMP-proton (WIMP-neutron) 

reduced mass, and Op™ are the proton and neutron 
cross section limits when a n , P — respectively. The cr| m 
is evaluated by attributing the entire counting rate to the 
nuclide A, so Op™ are overestimated by a factor f^ 1 , 
where f A is the isotopic fraction ( t " t ° al scnlmvcTncti ) 
of the isotope. Since YIa f A = 1> ana - f A x ; 
the constant of proportionality is taken as a more re- 
fined WIMP-proton (WIMP-neutron) cross section limit 

—lira 
p,n ' 



obtained from Eqn. (J2J with the relation 



Eqn. Q for the single proton (neutron) cross section 



Cp (cr„) reads: 
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Thus, given the cr p % ^ A \ model-independent limits on 
spin-dependent WIMP-nucleon interaction can be formu- 
lated either in terms of the a Pin or the coupling strengths 
a Pin . Although equivalent formulations, I discuss in a p , n , 
with both elaborated in a forthcoming paper. 

Substituting Eqns. £[J - into the obvious relation 
-t^- < 1 and taking m n m p , it can be shown that [5| 
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where the sign of the addition in parenthesis is that of 
the ^g"^ ratio. The allowed values of a p and a n are con- 
strained to the inside of a degenerate conic in the a p — a n 
plane, namely to a band between two parallel straight 
lines, whose stiffness is — $g"^ . This indicates the in- 
ability of a single nuclide experiment to fully constrain 
all tree level theoretical scenarios. 

For detectors with more than one active nuclide, as 
pointed out above, Op™^ are overestimated by a factor 
/J . Consequently the limits of Eqn. (@J become, for 
each nuclide, 
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Since J2a /a = 1 ) then: 
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where the sum is intended over the nuclear species 
present in the detector's sensitive volume. 

Eqn. © removes the conic degeneration of Eqn. Q, 
even if a nuclide that has little sensitivity to the WIMP- 
neutron interaction, e.g. Na or CI, is included. The 
allowed region is now a conic, and the exclusion results 
of multinuclide detectors are usually ellipses in the a p —a n 
plane. 

Since <r l p™ are functions of the WIMP mass M x 
only, Eqn. © is a three parameter relation, which yields 
a 3D allowed region, shaped as a tube with varying el- 
liptical cross section transverse to the mass axis. Results 
can be reported GHHSOl by plotting dp™^ vs M x , plus 
allowed regions for selected M x values. This is shown for 
several current experiments (Table 1) in Fig. ^ using 
spin values from Table 2. The calculation of Na and I 
for DAMA/Nal require recalculation of cr^ m for 
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Na and I from a mixed model |10| a p and cr„, followed 
by re-application of Eqns. J5J. While not completely 
rigorous, this is a good approximation. The DAMA/Xe- 
2 experiment, instead, uses almost pure 129 Xe, so it 
suffices to multiply the 0% reported in Ref. 
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CRESST did not 



report a an m< " A \ but since 16 is an even-even, spinless 
and doubly magic nucleus with no magnetic moment , 
its spin-dependent response has been neglected, and the 
exclusion of Ref. has been assumed a result of Al 
only. 



Experiments based on natural Ge and Si are not in- 
cluded, since their only spin-dependent sensitive isotopes 
are 73 Ge (7.8 % natural Ge) and 29 Si (4.67 % natural 
Si), respectively. Such small isotopic abundances lower 
the spin-dependent effective mass of detectors based on 
natural Ge and Si by a factor 10. 

The allowed a P; „ for the null result experiments of Fig. 
dare shown in Fig. 0for M x = 50 GeV/c 2 . This mass 
lies in the range 20 to 100 GeV/c 2 , where most experi- 
ments reach maximum sensitivity (smallest Op™^). 
Outside this range, the area of an experimentally allowed 
ellipse increases. Owing to the difference in sign between 
the ratios ^§3- of F and I, the angle between the ellipses 
of F-based and iodine-based experiments is large, allow- 
ing for an intersection significantly smaller than the orig- 
inal ellipses of each experiment. Although DAMA/Xe-2 
only provides two parallel lines, these are almost vertical 
due to the low value of < S p > , and the intersection of the 
DAMA/Xe-2 and the NAIAD ellipses is a very thin slice. 
Note that CRESST at 1.51 kgdy restricts a p at almost 
the same level as NAIAD, which has a much larger ex- 
posure (3879 kgdy). SIMPLE and PICASSO already ex- 
clude a small fraction of the NAIAD allowed ellipse, even 
with prototype low exposure results. Current limits from 
the intersection of all these experiments (shaded area) at 
M x = 50 GeV/c 2 are \a n \ < 1.3 and \a p \ < 1.5, and are 
effectively determined by only NAIAD and DAMA/Xe-2. 
In terms of cross sections, these are a n < 0.6 and a p < 0.8 
pb. Consequently, spin-dependent WIMP candidates in- 
teracting mainly with protons are excluded at about the 
same level of those interacting mainly with neutrons. As 
per the interaction strength, the limits on a p . n from null 
result experiments do not yet exclude a WIMP candidate 
with an ordinary weak interaction strength (|a p .„| = 1). 

If an experiment is equally sensitive, at a given M x , to 
all possible candidates, its ellipse becomes a circle. This 
can be seen by rewriting Eqn. © as: 
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where the coefficients a, f3, and 7 are given by: 

= J^A iim(A) 

P = Z A ± 
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The sign of each term in the /3-summation is given by the 
sign of the ratio ^g"^ for the isotope of mass number A. 
The coefficients a, (3, and 7 determine the shape of the 
ellipse in the a p — a n plane, hence its semiaxes and the 
angle between the major semiaxis and the a p axis. To 
make (3 = in Eqn. Q, corresponding to ellipse axes 
parallel to the coordinate axes, the detector must contain 
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FIG. 1: a) cr" m(A) and b) cr" m(A) vs M x for various elements by data from SIMPLE (thick solid), NAIAD (solid), CRESST 
(dotted), DAMA/Xe-2 (dash-dotted), DAMA/Nal (dotted for Na, dashed for I), and Tokyo (dashed). 



TABLE I: current experiment exposures (detector active mass times measurement duration). 



experiment 


NAIAD [1] 


DAMA/Xe-2 [2] 


DAMA/Nal [10] 


PICASSO [13] 


SIMPLE [4] 


Tokyo/NaF [3] 


CRESST [12] 


exposure (kgdy) 


3879 


1763.2 


57986 


0.056 


0.19 


3.38 


1.51 



TABLE II: Spin values for relevant nuclides. 



Nucleus 


Z 


J 


< S p > 


< S n > 


Ref. 


19 F 


9 


1/2 


0.441 


-0.109 


[5] 


23 Na 


11 


3/2 


0.248 


0.020 


[14] 


27 Al 


13 


5/2 


0.343 


0.030 


[15] 


35 CI 


17 


3/2 


-0.059 


0.011 


U 


37 C1 


17 


3/2 


-0.178 


" 


[4] 


127j 


53 


5/2 


0.309 


0.075 


[14] 


129 Xe 


54 


1/2 


0.028 


0.359 


[14] 



"cal cula ted in the odd group approximation |4|, using data from 
Ref. Il J . Since 37 CI and 39 K have the same number of neutrons 
and similar number of protons, 37 C1 spin values can be evaluated 
more accurately than in Ref. Q by assuming the same neutron 
group spin and angular momentum as 39 K. 

nuclei with different signs of -^§3-; to get a = 7 both 
odd Z and odd N nuclei are needed. 

Until a positive signal is found, experimental improve- 
ments simply shrink the various ellipses. In the case of 
a positive signal, like DAMA/Nal, both an upper (a 1 ™ 1 ) 
and a lower (cr^ mm ^) limit on a a exist: 
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FIG. 2: Excluded regions in the plane for M x = 

50 GeV/c 2 , for the experiments of Table The resulting 
allowed region is determined by the intersection of NAIAD 
and DAMA/Xe-2. 



< f A < (9) 
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where the isotopic fractions /a are again due to the over- 
estimate of attributing the entire counting rate to iso- 



topic specie A only, which affects both a a" 1 and <j^ mm ^ 
by a factor f^. Analogy with Eqn. @ then leads to the 
auxiliary cross sections: 
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FIG. 3: Permitted regions for a positive signal (elliptical 
shell), and two null signal (dashed and solid) experiments. 
The dotted inner lines are arbitrary and for illustrative pur- 
pose only. The 2 thick border areas are the intersection of cur- 
rent DAMA/Nal and DAMA/Xe-2 results; if also DAMA/Xe- 
2 had a positive signal this intersection would become the 4 
hatched regions; a positive signal F-based experiment could 
reduce such intersection to the two solid regions (see text). 
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allowing to complete Eqn. @ with: 
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These extra equations transform the a p — a n plane ellipses 
into elliptical " shells" and pairs of bands symmetric with 
respect to the origin. WIMP masses are excluded when 
the lower limits of Eqn. are incompatible with the 
upper limits of Eqns. 

Fig. |2| shows the intersection (thick border) of 
DAMA/Nal and DAMA/Xe-2 permitted regions for 
M x = 50 GeV/c 2 . NAIAD is reaching this region but 
to date cannot cut any part of it. When the positive 
signal of DAMA/Nal is considered, the limits for M x = 
50 GeV/c 2 become < \a n \ < 1.3, 0.5 < \a p \ < 1.1, or 
< <r„ < 0.6, 0.07 < cr p < 0.4 pb. Though tighter than 
the above, these limits on the WIMP-nucleon coupling 
strengths are still compatible with the ordinary weak in- 
teraction strength. 

The two allowed areas are symmetric with respect to 
the origin and correspond to one in the cross section 
representation. This symmetry and doubling are due to 



the fact that a Pi „ cx ±y/Op~^, so that each point in the 
cross section representation corresponds to two symmet- 
ric points in the coupling strength representation. 

The inability to set a lower limit for \a n \ is due to 
the null signal of the DAMA/Xe-2 experiment. Had 
DAMA/Xe-2 a positive signal, its permitted region would 
become a pair of bands symmetric with respect to the ori- 
gin of the a p — a n plane, leading to the four intersection 
regions (two in the cross section representation) shown 
in Fig. |3 as hatched. To reduce the number of allowed 
regions, a third experiment with a different orientation is 
needed, such as one of the F-based experiments of Fig. [2] 
Fig. |21 also contains a projection for PICASSO (dashed), 
along with an hypothetical positive signal (dotted) ; other 
F-based experiments require improvements or higher ex- 
posures to reach the same level. Note that F-based ex- 
periments seem competitive even with exposures signifi- 
cantly lower than sodium iodide-based. 

In the case of Fig. [21 the F-based experiments would 
remove completely two of the hatched areas, allowing 
only the two solid regions symmetric with respect to the 
origin. The two regions correspond to positive (repulsive) 
and to negative (attractive) WIMP-nucleon interaction. 

If there is a single WIMP specie, there must be a pair 
of symmetric regions (one in the cross section representa- 
tion) permitted by all experiments with positive signal. If 
there is no such common permitted region, the existence 
of a single WIMP specie is excluded. 
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